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In this paper some specific issues related to point defects in GaAs nanowires are addressed witfi 
tfie fielp of density functional theory calculations. These issues mainly arise from the growth of 
nanowires under conditions different from those used for thin films or bulk GaAs, such as the 
co-existence of zincblende and wurtzite polytypes, the use of gold particles as catalyst, and the 
arsenic-limited growth regime. Hence, we carry out density-functional calculations for As vacancies, 
GaAs antisites, and Au impurities in ZB and WZ GaAs crystals. Our results show that As vacancies 
can diffuse within in a ZB GaAs crystal with migration barriers of ^1.9 eV. Within WZ GaAs, As 
vacancy diffusion is found to be anisotropic, with low barriers of 1.60 up to 1.79 eV (depending on 
doping conditions) in the a6-plane, while there are higher barriers of 2.07 to 2.44 eV to diffuse along 
the c-axis. The formation energy of Au impurities is found to be generally much lower than those 
of arsenic vacancies or GaAs antisites. Thus, Au impurities will be the dominant defects formed 
in Au-catalyzed nanowire growth. Moreover, we find that it is energetically more favorable by 1 
to 2 eV for an Au impurity to replace a lattice Ga atom than a lattice As atom in GaAs. An Au 
substitutional defect for a lattice Ga atom in ZB GaAs is found to create a charge transfer level in 
the lower half of the band gap. While our calculations locate this level at Ey + 0.22 eV, taking into 
account the inaccuracy of the density functional that ought to be corrected by a downshift of Ey by 
about 0.2 eV results in good agreement with the experimental result of Ey + 0.4 eV. 
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I. INTRODUCTION 

Semiconductor nanowires have been emerging as a 
promising building block for various nano-devicesA These 
potential applications include light emitters P solar cells'^ 
and microelectronics.^ ^ All these applications require 
such a prerequisite that it is possible to keep the con- 
centration of intrinsic defects and impurities in the 
nanowires below a given threshold so as to enable long 
exciton lifetimes and to avoid unintentional doping due 
to electrically active impurities. In recent years, many 
experimental efforts^ have been devoted to fast and 
controllable growth of nanowires that is desired for re- 
liable devices; and a gold nanoparticle is often used as 
a catalyst to promote the growth of GaAs nanowires. 
However, an Au droplet that accelerates an underneath 
growing nanow ire ca n also leave impurities within the 
GaAs nanowire.H^ Previous experimental efforts have 
focused on characterizing Au defect levels within a GaAs 
crystal Moreover, experimental studies employing 
X-ray energy-dispersive spectroscopy^ indicate devia- 
tions from stoichiometry in the nanowires close to the 
growth zone, and thus point to the abundant intrinsic 
defects, in addition to Au impurities, in GaAs nanowires. 
While carrier lifetimes in core- shell nanowires from Au- 
free self-catalyzed growthF^^ were found to be much 
longer than in nanowires from Au-catalyzed growth^^, 
they still fall short of the values reported for bulk sam- 
ples by several orders of magnitude. This may indicate 
a higher level of growth-related intrinsic defects (as com- 
pared to bulk) even in Au-free nanowires. 

The focus of our interest is the defects that could be 



formed in nanowire growth below the metal nanoparti- 
cle. While the most abundant defect in low-temperature 
grown bulk GaAs is the As antisite^^, we expect this de- 
fect not to play a role under the more Ga-rich growth 
conditions beneath a metal catalyst particle. It is known 
that near the melting point and under the conditions typ- 
ical of liquid-phase epitaxy. As vacancies are the domi- 
nant defects in GaAs.^^ Furthermore, under As-deficient 
conditions, the As lattice site might remain unoccupied 
(As vacancy Vas), it could be occupied by gallium (Ga 
antisite GaAs), or, in case of growth with a gold catalyst, 
may be occupied by an Au atom (Auas)- Since there is 
a general interest in possibly detrimental effects of Au 
on the properties of GaAs nanowires, we include AuQa 
defects into our study. 

In this context, one has to consider that the effective 
growth conditions for the nanowire material are probably 
very different from the moderately arsenic-rich conditions 
commonly used in GaAs thin film growth. This is due 
to the following reasons: In nanowire growth, material 
deposited from the vapor phase onto the substrate, the 
sidewalls and the metal nanoparticle can reach the in- 
terfacial area between the nanowire tip and the particle 
via surface diffusion^ 20 22 diffusion through the liq- 
uid metal particle.'^ It is clear that sufficient Ga atoms 
can reach the interfacial growth zone at the nanowire 
tip, since the catalyst particle actually consists of a Au- 
Ga alloy. ^ The situation is less clear concerning the As 
supply, since the low-pressure solubility of As in gold is 
low (in Au-catalyzed growth), or a high high tempera- 
ture leads to As loss (in self-catalyzed growth). In these 
cases, one would expect As-deficient growth conditions 



at the nanowire-particle interface. One possible scenario 
is that additional arsenic can reach the growth zone via 
an As vacancy diffusion mechanism through the GaAs 
nanowire. Moreover, under As-deficient conditions, one 
would also expect that As vacancies could be left be- 
hind within the GaAs nanowire by the advancing inter- 
facial growth zone. A subsequent annealing process in 
arsenic vapor might be helpful to remove these vacan- 
cies. Therefore, it is important to understand Vas diffu- 
sion within GaAs bulk, in order to comprehensively un- 
derstand the As supply path during the growth process 
and to estimate a suitable annealing temperature for re- 
moving VasS in a nanowire. A previous t heoret ical study 
based on density functional theory (DFT)P^^ calculated 
relatively high migration barriers of ~ 2.4 eV for a Vas in 
zincblende GaAs.^^ This seemingly implies that As trans- 
port through the nanowire could be inefficient. However, 
the convergence of the atomic-orbital basis sets employed 
in these DFT calculations has not been demonstrated. 

While wurtzite GaAs cannot be obtained in bulk form 
by ordinary growth techniques, the GaAs nanowires 
grown with Au catalyst may exhibit both zincblende 
(ZB) and wurtzite (WZ) structures, or alternating seg- 
ments of these (and other) polytypes through stacking 
fault s.f^l^^llll Hence, it is crucial to understand Vas, GaAs, 
and Au defects in both ZB and WZ GaAs crystals. There 
is consider able k nowledge, both from experimental^^ and 
theoreticaP^I^ sides, about intrinsic point defects in ZB 
GaAs. However, experimental data on point defects in 
WZ GaAs is still elusive. In this work, we will perform 
plane- wave DFT calculations to study various properties 
associated with Vas, GaAs, and Au defects in ZB and WZ 
GaAs crystals, and to investigate and compare Vas dif- 
fusion in both GaAs polytypes. This paper is organized 
as follows: The computational approach and supercell 
models are described in Sec. |n| The formation energies 
of Vas, GaAs defects and Au impurities in ZB and WZ 
GaAs crystals are discussed in Sec. 
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VasS in various charge states in ZB and WZ GaAs crys- 



tals are presented in Sec. IV We summarize our work in 
Sec. El 



II. METHODS OF CALCULATION 

In this work, we have performed DFT calculations to 
study Vas, GaAs and substitutional Au impurities in ZB 
and WZ GaAs using the supercell method. The con- 
ventional unit cell for a ZB and a WZ GaAs crystal 
is shown in Fig. [T] Sufficiently large supercells of 216 
atoms (ZB) or 96 atoms (WZ) are used to fully relax the 
strain induced by defects and impurities. We employ the 
projector augmented- wave (PAW) method^^ as imple- 
mented in the VASP codd^^ for all DFT calculations. 
The PAW potentials include 454p, 4^4^, and 5d6s elec- 
trons as valence electrons for As, Ga, and Au, respec- 
tively. The generalized-gradient approximation (GGA) is 
used for the exchange-correlation functional. The Bril- 




FIG. 1: (Color online) The structure of ZB and WZ GaAs 
crystals. The conventional unit cell for a ZB and a WZ GaAs 
crystal is shown in (a) and (b), respectively. The As and Ga 
atoms are represented by large pink and small green spheres, 
respectively. 



louin zone integration is performed using the Monkhorst- 
PaclP^ scheme over a (2 x 2 x 2) k-point mesh (shifted 
with respect to the F-point) for both ZB and WZ su- 
percells. Then, the partition length is 0.06 Bohr~^ or 
smaller in each direction. Spin-polarized calculations 
are performed for supercells that contain an odd num- 
ber of electrons. The optimized kinetic energy cutoff of 
^cut = 250 eV is used to compute the formation energies 
of various impurities and defects, whereas an Ecut of 200 
eV is used for studying Vas diffusion in GaAs. The op- 
timized parameters allow us to obtain a relative energy 
convergence of within 10 meV for all systems presented 
in this study. The method of a homogeneous background 
charge is used to model charge states of the defect within 
the supercell approach. We determined various diffusion 
paths for VasS within the As sublattice of both ZB and 
WZ GaAs crystals, employing the nudged-elastic band 
(NEB) method. ^^ ^^ Here, the migration pathway is rep- 
resented by several intermediate configurations between 
two fully relaxed end points and each image is relaxed 
until the perpendicular forces with respect to the mini- 
mum energy path are less than a given tolerance, which 
is set to be 0.03 eV/A in our calculations. 

For ZB GaAs (space group F43m), we have obtained 
the optimized lattice constant of a = 5.76 A, in good 
agreement with the experimental value of a = 5.653 A.l^ 
For WZ GaAs (space group P63mc), we have computed 
the lattice constant to be a = 4.05 A and c = 6.69 
A, which compares well with previous GGA results of 
a = 4.050 A and c = 6.678 AF^ The calculated direct 
band gaps of GaAs are 0.17 eV and 0.22 eV in ZB 
and WZ crystals, respectively. These results refiect the 
well-known underestimation of the electronic band gap 
by semi-local DFT functionals, in particular, when the 
(too large) theoretical lattice constant is used. In addi- 
tion, we have performed DFT calculations with the same 
geometries as in GGA, but using the hybrid functional 
proposed by Heyd, Scuseria, and Ernzerhof (HSE).S^Our 
calculations include 25% of the exact exchange, and the 
screening parameter is set to be 0.2 A~^. The band gaps. 
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that are computed to be 1.15 eV and 0.99 eV in ZB and 
WZ GaAs, respectively, are much closer to experimental 
values. Moreover, we conclude from these calculations 
that the valence band maximum experiences a down-shift 
of 0.48 eV in ZB GaAs (0.36 eV in WZ GaAs). These val- 
ues are in good agreement with recent literaturd^, and 
will be used as corrections when our computed results are 
compared to experiment. The orthorhombic Ga (space 
group Cmca) crystal and trigonal As (space group R3m) 
crystal are used to compute the Ga and As chemical po- 
tentials. 

Since ZB and WZ segments co-exist in GaAs 
nanowires, the band edges from ZB and WZ supercell 
calculations should be aligned accordingly. We have con- 
structed a 50 atom supercell of a heterostructure of WZ 
and ZB GaAs. The supercell consists of 13 ZB and 12 WZ 
bilayers that stack along [111] or [0001] for a ZB or a WZ 
segment, respectively. The constructed heterostructure 
is optimized for both the lattice constant and internal 
coordinates. Note that the lateral dimension of the het- 
erostructure happens to coincide with the average of the 
WZ and ZB lattice parameters. The valence band offset 
of the ZB and WZ segments can be computed^^ as 



AEvsiGsiAs) = £ 



ZB 



+ (^f<f,r*5(As) - EfJ, -s(As)), (1) 

where E'^^ and E'^'^ are the valence band top en- 
ergies at the F-point for ZB and WZ GaAs, respec- 
tively, £^^^g(As) and E^^^{As) are the core level ener- 
gies of As atoms in bulk ZB and WZ GaAs, respectively, 
^co?e'^^(As) and E^^/'^^iAs) are the core level energies 
of As atoms in ZB and WZ segments of the heterostruc- 
ture, respectively. We compute the band offset to be 
AEvsiGdiAs) = —0.0601 eV, in good agreement with 
the value of —0.0632 eV from a previous calculation.!^ 




-0.4 -0.2 0.0 0.2 



Fermi level - 




Fermi level 



(eV) 



FIG. 2: (Color online) The formation energy of an As vacancy 
(a) and a Ga antisite (b) in various charge states in GaAs 
polytypes at different doping conditions, as a function of the 
Fermi level sf- The formation energies of As vacancies in 
ZB GaAs and WZ GaAs are shown as dashed and solid lines, 
respectively. The thick lines represent stable charge states. 
The vertical lines indicate the HOMO and LUMO levels in 
the defect-free supercell and thus delimit the energy interval 
where the GGA total energies are considered reliable. The 
left and right ^-axes show the formation energies under As- 
rich (/i^g^^^) and Ga-rich (/Xa^^^) conditions as described in 
the text. Here, one has the identity /i^^,^^^ = /i^f """^ - 0.7 eV. 



III. THE ENERGETICS OF DEFECTS AND Au 
IMPURITIES IN GaAs NANOWIRES 

We start by calculating the formation energies of rel- 
evant intrinsic defects. The formation energy of a q- 
charged {q is an integer number) defect (or impurity) 
X at site F as a function of the Fermi level is defined 

by 

Ef [X^,eF] = E [Xl.]-Ehost-J2ni^ii+q{eF+E,)+E,„rr 

i 

(2) 

Here, E [Xy] is the total energy of a g-charged supercell 
with a defect Xy and Ehost is the energy of a perfect 
GaAs bulk supercell with the same shape, rii and fii 
represent the change in the number of species i in the 
supercell and the chemical potential for a species z, re- 
spectively, is the Fermi level that depends on the 
doping condition, and Ey is the valence band top en- 
ergy of the corresponding bulk GaAs crystal. Due to the 
(2x2x2) k-point mesh used for supercell calculations, the 



gap between the highest occupied (HOMO) and the low- 
est unoccupied (LUMO) orbitals is much higher in these 
supercell calculations compared to the gap at the F-point 
in GaAs bulk calculations. Even though this widened gap 
is a technical artifact of our method, it puts us in position 
to vary ep in a rather wide range while still maintaining 
the physically correct charge distribution localized near 
the defect. The energy interval around Ey in which our 
calculations can be expected to yield physically mean- 
ingful charge transfer levels, as a charged defect cannot 
transfer charge to the LUMO or from the HOMO, is indi- 
cated in the figures. The linear and quadratic correction 
term with respect to q of Ecorr = AVq + aq^ /2eL ac- 
counts for the spurious electrostat ic int eraction between 
periodic images of charged defect AU is a correc- 
tion for the spurious potential off-set induced by the finite 
defect concentration in the calculation; it is calculated 
from the energy difference of atomic core levels between 
a neutral supercell with an Xy defect and a perfect bulk 
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supercell. The values for a and e, the Madelung constant 
and the static dielectric constant, respectively, are taken 
according to the ZB or WZ polytype of bulk GaAs, and 
L is the supercell dimension used in the corresponding 
calculations. The values of a, e, and L are listed in Ta- 



ble III in the Appendix for the cases of ZB and WZ GaAs 
crystals, respectively. Our computed dielectric constant 
e for ZB GaAs is 32% larger than the experimental value 
of 0.90 eA~^V~^.^^ At the growth condition, the chemi- 
cal potentials of Ga and As are in equilibrium with bulk 
GaAs. Thus, one has the identity /iQa + A^As = MCaAs^ 
where /icaAs = ~^GaAs cohesive energy per for- 

mula unit of a pertinent GaAs polytype. The range of 
the As chemical potential is specified in accordance with 
the growth condition, with the upper bound set to be at 
the equilibrium with bulk As, i.e., /i^f = -E^j^J^. The 
lower bound is set to be at the equilibrium with bulk Ga 
as /i^,^^^ = -^GaAs + ^Ga^- Disregarding the difference 
of the cohesive energies of GaAs of 20 meV in ZB and 
WZ crystals, we align the upper bound ja^f^^ for ZB 
and WZ to the same value. We take the lower bound 
flower to be 0.7 eV below /i^f so that /i^^^^^ is very 
close to a Ga-rich condition for both ZB and WZ GaAs. 
The bulk energies of GaAs, Ga and As crystals are calcu- 
lated using the optimized lattice parameters. Note that 
the upper and lower bounds of the As chemical poten- 
tials correspond to As- and Ga-rich growth conditions, 
respectively. In addition, the cohesive energy bulk Au 
per atom is inserted for the Au chemical potential, i.e.. 

Figure |2|a) shows the Vas formation energies for var- 
ious charge states indicating the formation energies 
of Vas in WZ are overall lower than those in ZB. The 
formation energies are computed using Eq. [2j and the 
potential off-sets AV associated with the VasS are listed 
in Table |IV| in the Appendix. Under an As-rich growth 
condition (left ?/-axis in Fig. [2|a)), the formation ener- 
gies of neutral VasS are quite high (3.3 eV and 3.2 eV 
in ZB and WZ, respectively). Under a Ga-rich condi- 
tion (right y-ax.is in Fig. [2ja)), the formation energies 
are lowered by 0.7 eV; however, the energies remain rel- 
atively high. Thus, the formation of As vacancies could 
be energetically unfavorable under usual GaAs nanowire 
growth conditions in both ZB and WZ segments. Our re- 
sults show that and V^^ defects are stable at various 
doping conditions in ZB and WZ GaAs crystals. Deep 
(+/-) levels are found at +0.44 and +0.45 eV 
in ZB and WZ, respectively. Since GaAs nanowires may 
have coexisting ZB and WZ segments, it is instructive 
to specify the charge transfer levels of defects in the WZ 
segments also with respect to the valence band edge of 
the ZB segments, using our results that the valence band 
maximum of ZB GaAs is 0.0601 eV lower than in WZ 
GaAs. Thus, one obtains the (+/— ) level in WZ to be 
E^^ + 0.51 eV. We found that V^^ and V^^ in both 
ZB and WZ GaAs crystals are unstable, consistent with 
Ref. 30, because the structures relax to a GaAs antisite 
that is neighboring a Ga vacancy. As seen from Fig. |2j 



the neutral As vacancies are energetically unfavorable, 
indicating that a Vas in GaAs is a so-called 'negative 
/7' system. Our results are consistent with a previous 
theoretical studie d^^ ' ^^'^ -*^, although the quoted absolute 
positions of the (+/— ) levels in these studies are different 
from ours. For a quantitative comparison, it is important 
to realize that Schulz and von LilienfelcF^ attempted to 
set Ey to the true ionization potential, whereas we use the 
plain DFT-GGA value. Estimating the true valence band 
top from our HSE calculations, a down-shift of 0.4 eV 
(see Section |ll|) should be applied to our Ey value in or- 
der to compare to their work, or to experimental data. 
The same applies when comparing to the values of Komsa 
and Pasquarello,'^ who performed all their calculations 
with the HSE functional. With this in mind, our results 
are in much better agreement with these previous works. 
It is worth noting that the negative U feature of a Vas 
disappears in their HSE calculations. This is due to the 
fact that the (negative) contribution of the electronic ex- 
change energy, which is more pronounced in the hybrid 
functional, energetically favors spin-polarized solutions, 
e.g., the neutral Vas in the present case. Whether this 
stabilization of a V^g in a narrow range of the Fermi level 
is indeed a physical feature or not should be determined 
by future experiments. 

The Ga antisite GaAs could be an another intrinsic de- 
fect within a growing GaAs nanowire under As-deficient 
conditions, and we show the formation energies of GaAs 
as a function of in Fig. [2|b). The formation ener- 
gies are computed using Eq. [2] with the potential off-set 
A V associated with GaAs defects listed in Table [IV| in 
the Appendix. The formation energies are specified for 
both Ga-rich (right y-ax.is in Fig. [2|b)) and As-rich (left 
^-axis Fig. [2|b)) conditions, and formation energies are 
1.4 eV lower in a Ga-rich condition than in an As-rich 
condition. In both ZB and WZ GaAs crystals, GaAsS 
with charge states from 2+ to 2— may be stable under 
certain doping conditions. The associated formation en- 
ergies in WZ GaAs are overall lower than in ZB GaAs. 
The Ga^g and GaJ^ defects are energetically more favor- 
able by 0.13-0.14 eV in WZ than in ZB, whereas a Ga^^ 
defect is only 0.01 eV more stable in WZ. As a result, 
the charge transfer levels (0/—) differ by 0.12 eV in ZB 
and WZ. It follows that, on an absolute scale, the charge 
transfer levels are located at £;^^+0.21 and 0.33 eV in ZB 
and WZ segments of a GaAs nanowire, respectively. We 
notice that both intrinsic defects of Vas and GaAs have 
deep charge transfer levels and they are deeper in WZ 
GaAs. Comparing Fig. [2|a) and (b), one may find that 
VasS are more stable than GaAsS for the doping condi- 
tions of Sf - E^^ < 0.21 (0.27) eV in ZB (WZ) under As 
rich conditions, while GaAsS are more stable than VasS for 
the doping conditions of ep-^^^ > 0.21 (0.27) eV in ZB 
(WZ). When comparing to experiment, one should take 
into account that the "true" valence band top is lower 
than the Ey resulting from our DFT-GGA calculations, 
such that the range of horizontal-axis values plotted in 
Fig. [2] approximately reflects the "true" band gap. 
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FIG. 3: (Color online) The formation energies of a substitu- 
tional Au defect in the ±l-charged, it2-charged, and neutral 
states in WZ and ZB GaAs crystals. The formation ener- 
gies of a substitutional Au impurity for a lattice As atom 
and a lattice Ga atom are shown in (a) and (b), respectively. 
The formation energies of Au impurities in ZB GaAs and WZ 
GaAs are shown as dashed and solid lines, respectively. The 
thick lines represent stable charge states. The vertical lines 
indicate the HOMO and LUMO levels in the defect-free su- 
percell, same as Fig. [2] The left and right ^-axes show the for- 
mation energies under As-rich (a^as^^^) Ga-rich (/Xa^^^^ 
conditions as described in the text. 



We also study the energetics of substitutional Au im- 
purities in GaAs crystals, and the computed formation 
energies of Auas and AuQa as a function of ep are shown 
in Fig. |3] The formation energies of Auas and AuQa 
are computed using Eq. |2j with the potential off-sets as- 
sociated with Auas and AuQa listed in Table |IV| in the 
Appendix. Moreover, the formation energies are speci- 
fied for both Ga-rich (right y-ax.is in Fig. |3| and As-rich 
(left 7/- axis Fig. [3| conditions. It is found that under an 
As-rich growth condition an AuQa is energetically more 
favorable than an Auas by about one to two eV both 
in ZB GaAs and in WZ GaAs, depending on the dop- 
ing conditions. However, under a Ga-rich growth condi- 
tion, the formation energy of an AuQa is 0.7 eV higher 
than that under an As-rich condition, while the forma- 
tion energy of an Auas becomes 0.7 eV lower. Therefore, 
under a Ga-rich condition, the Auas defects turn out to 
be more stable than AuQaS. We expect that Auca de- 
fects with charge states from 2+ to 2— may be stable 
under certain doping conditions. An AuQa has its charge 



transfer levels mostly in the lower part of the band gap, 
while the charge transfer levels of an Auas lie at some- 
what higher energies. The Auas impurities occur pref- 
erentially as 2+, -f or neutral defects, while the levels 
of negatively charged defects lie probably already above 
the conduction band minimum of the host material. Fig- 
ure |3|a) indicates the charge transfer levels of AuasS are 
deep, E^^ + 0.78 eV and E^^ + 0.91 eV for ZB and 
WZ crystals, respectively. In contrast, an AuQa impu- 
rity switches from a neutral state to a 1— charged state 
at +0.22 eV and +0.35 eV in ZB and WZ crys- 
tals, respectively, as seen in Fig. [sj^b). As stated before, 
one should allow for a correction of about 0.4 eV to our 
Ey value when comparing to experiment. Experimen- 
tal studies using deep-level transient spectroscopy and 
photoluminescence spectroscop}'^^ have identified an 
Au-related deep acceptor level of about 0.4 eV above the 
valence band in bulk GaAs. Moreover, we find an AuQa 
has a lower formation energy compared to an Auas under 
the moderately As-rich growth conditions conventionally 
used. Taking into account the down-shift of Ey with re- 
spect to the DFT-GGA value (see e.g., Ref. |44|, it is 
plausible that the experimentally observed defect level 
at 0.4 eV above the valence band is indeed due to AuQa 
defects. 



We conclude this Section by discussing the energetics 
of defects caused by a deficiency of arsenic under varying 
growth conditions. In Fig. [4]^ a), we compare the forma- 
tion energies of the lowest energy Vas, Auas and GaAs 
defects in p-type ZB and WZ GaAs crystals with the 
condition of ep — E^^ = 0, as a function of As chemical 
potential /iAs- Our calculations show that Auas impu- 
rities have the lowest formation energies, while Vas and 
GaAs defects have clearly higher formation energies. Fig- 
ure [Ij^b) shows the same comparison for n-type ZB and 
WZ GaAs crystals with = + 1.0 eV. Here, a GaAs 
in the 2— charge state has the lowest formation energy, 
followed by Auas and Vas- Consequently, when grow- 
ing GaAs nanowires with a gold catalyst droplet, some 
Auas defects will be formed in the GaAs nanowire, in 
particular for a p-doped material. The Au impurities 
act as deep centers and are thus detrimental to the op- 
tical properties of nanowires. In gold-free self-catalyzed 
growth, GaAs and Vas defects may exist in the nanowire, 
albeit at much smaller concentration than the Au de- 
fects if the material is p-doped. Out of these three de- 
fect species, only the Vas defect is expected to be mobile 
via a hopping diffusion mechanism. The mobility of the 
GaAs and Auas defects will most probably be vacancy- 
mediated, and thus require the Vas presence. Moreover, 
it is conceivable that Vas hopping diffusion could play a 
role in the nanowire growth process as a mechanism sup- 
plying arsenic to the interfacial growth zone between the 
nanowire tip and the catalyst droplet. Hence, we will in- 
vestigate the detailed diffusion processes of VasS in GaAs 
crystals in the following Section. 
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FIG. 4: (Color online) The formation energies of Vas, Auas 
and GaAs defects in GaAs crystals as a function of /xas- In 
(a) and (b), the Fermi energy is set to be at the calculated 
top of valence band and 1 eV above, respectively, to describe 
a p-type or an n-type material. Only the energetically most 
favorable charge states are shown. 



IV. As VACANCY DIFFUSION IN GaAs 




FIG. 5: (Color online) The geometry of VasS in GaAs crystals. 
The relaxations of and in ZB GaAs are illustrated in 
(a), while the relaxations of two sets of stable and in 
WZ GaAs are characterized in (b) and (c), respectively. We 
have found two stable V^^s (Ml and M2) and two V^^s (Ml 
and M2) in WZ GaAs. A Vas is highlighted by a big blue 
sphere, and its four neighboring Ga atoms are highlighted by 
A, B, C, and D labeled spheres. These neighboring Ga atoms 
form a tetrahedron. The other attached As and Ga atoms 
are represented by small pink red and large green spheres, 
respectively. The right panel indicates the relaxed edge length 
of the Ga tetrahedron. The upper and lower numbers refer to 
a and a V^^, respectively. 



Arsenic vacancies in GaAs have been characterized ex- 
perimentally by positron annihilation^^ and by scanning- 
tunneling microscopy^. For ZB GaAs, a number of the- 
oretical st udies, using (semi-)local DFT^^ or hybrid 
functionalPII^, have been carried out, and the diffusion 
pathway of the arsenic vacancy has been investigated^^. 
Interestingly, these calculations claimed the existence of 
a metastable interstitial state along the diffusion path in 
ZB GaAs, while this issue is unexplored for WZ GaAs. 

We start by characterizing the geometry of VasS in 
GaAs crystals in more detail. As shown in Fig. [5j the 
neighboring Ga atoms with respect to a Vas define a Ga 
tetrahedron. To specify the local deformation of a crystal 
due to the V^s presence, we study the relaxation of the 
edge length / of the tetrahedron and compare it to the 
equilibrium distance between two neighboring Ga atoms 
in the ZB (WZ) GaAs crystal of /q = 4.07 (4.05) A. Simi- 
larly, the relaxed tetrahedron volume V will be compared 







yL 


^a"s 


ZB 


0.9221^0 


0.757^0 


0.516% 


WZ-Ml 


0.896^0 


0.696% 


0.518% 


WZ-M2 


0.926Vb 




0.512% 



TABLE I: Ga tetrahedron volumes of %vsS in various charge 
states. % represents the equilibrium tetrahedron volume in 
a ZB (WZ) GaAs crystal of % = 7.97(7.93) 



to the equilibrium tetrahedron volume in ZB (WZ) GaAs 
crystal of Vb = 7.97 (7.93) A^ 

For a in ZB GaAs, the edge lengths of the Ga- 
atom tetrahedron are all contracted (see the right panel 
of Fig. [sja)), resulting in a significantly reduced volume 
of F = 0.516%) (48.4% compression) as listed in Table [l| 
Especially, the A-D and B-C edges are most contracted 
up to / = 0.702/0 (29.8 % compression), while other edges 
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are compressed by ~11 %. The deformation reduces the 
symmetry of the distorted Ga tetrahe dron to D2d which 
is in accordance with previous studies.l^^In WZ GaAs, 
there exist two different, but energeticahy very similar 
structures (Ml and M2) for a relaxed V^^ (see Fig. [sj^b) 
and (c)). The WZ-M2 configuration is slightly more sta- 
ble than the WZ-Ml configuration in terms of energet- 
ics. Both WZ-Ml and WZ-M2 F^s configurations exhibit 
similar deformations as a V^^ in ZB, and the relaxed 
tetrahedron volumes are reduced by 48.4 % to 48.8 %, 
respectively. This is comparable to the ZB case. The re- 
laxed Ga tetrahedron for WZ-Ml and WZ-M2 V^^s pos- 
sesses a C2v and a Ci symmetry, respectively. 

For a the tetrahedron of neighboring Ga atoms 
contracts in volume by only 7.8 % in ZB GaAs, whereas 
it contracts by 7.4 % to 11.4 % in WZ GaAs, as listed 
in Table [ij In particular, V^^s have one edge compressed, 
while the other edges are slightly expanded. The WZ- 
M2 configuration has the A-C edge (within the ab 
plane) compressed. In contrast, the WZ-Ml V^^ con- 
figuration has the A-D edge (out of the ab plane) com- 
pressed. In other words, the deformation of in WZ 
GaAs is anisotropic, and hence the relaxed Ga tetrahe- 
dron for WZ-Ml and WZ-M2 V^^s possesses C2v or Ci 
symmetry, respectively. In addition, despite the neutral 
state being unstable, we mention that the volumes of 
neutral As vacancies are 0.757Vo in ZB and 0.696Fo in 
WZ-Ml, respectively. It is worth noting that the va- 
cancy volume becomes smaller with increasing electronic 
charge. The physical reason for this general trend lies in 
the fact that the additional electrons in the vacancy al- 
low for the filling of orbitals with bonding character that 
are formed by symmetry-adapted linear combinations of 
the Ga dangling-bond orbitals. 

In the following, we assume that the Vas diffusion pro- 
ceeds via a neighboring As atom of the arsenic sublattice 
hopping into the vacancy. Schematically, an As atom 
hops out of an initial Ga-tetrahedron cage into an adja- 
cent target Ga-tetrahedron cage that previously enclosed 
the As vacancy site. We refer to the tetrahedron that sur- 
rounds the migrating As atom as the initial tetrahedron, 
and refer to the adjacent Ga-tetrahedron that surrounds 
the target As vacancy site as the target tetrahedron. We 
note that the reverse process occurs with the same prob- 
ability, with the role of initial and target tetrahedron 
interchanged. The diffusion process can be decomposed 
into an As escape from the initial tetrahedron, its inter- 
stitial motion between the two tetrahedra, and an entry 
into the target tetrahedron. Note that one of the Ga 
atoms is located at the common apex of both tetrahedra. 

Since ZB GaAs is a cubic crystal, the four possible dif- 
fusion paths of a Vas to four neighboring As sites are sym- 
metrically equivalent. The diffusion path shown Fig. |6] 
is sufficient to describe As diffusion in ZB GaAs. Fig- 
urepfa) illustrates the detailed diffusion process of a 
in ZB GaAs, indicating an overall migration barrier of 
1.93 eV. The local minimum configurations {i and v in 
Fig. [oj^a)) as well as the saddle point configurations {ii 



(a) Vl/ZB 

2.5 l ■ ■ ■ ■ 




FIG. 6: (Color online) The energy path of V^^ (a) and V^' 
(b) diffusion in ZB GaAs. The structures of two ground state 
minima (i and v), an interstitial state (in), and transition 
states {ii and iv) are illustrated. The As and Ga atoms are 
represented by lager pink and small green spheres, respec- 
tively. 



and iv in Fig. [oj^a)) along the diffusion path are illus- 
trated by the ball-and-stick models in the figure. In the 
target tetrahedron, the AG edge corresponds to the short 
edge (the AD in Fig. |5|a) ). By expanding the bond be- 
tween As and Ga at the D-site, the diffusing As atom 
escapes from the initial tetrahedron. During the escape 
process, the initial tetrahedron is expanding, while the 
target tetrahedron is contracting. After passing through 
the ABC face of the tetrahedron, the diffusing As atom 
arrives at the first saddle point {ii) which represents the 
transition state for migration. In this configuration, 
the As-Ga bonds to the initial tetrahedron apexes are 
broken except for the A-site. The target tetrahedron 
volume is compressed from 0.922Vb to 0.719Fo which cor- 
responds to the volume of a neutral As vacancy (see Ta- 
ble [l|). This observation implies that a positively charged 



As vacancy in the target tetrahedron attracts electrons 
ahead of the As migration. When migrating toward the 
target tetrahedron in the interstitial space, the diffusing 
As atom reaches an interstitial state {Hi) with a shallow 
energy dip. At this moment, the diffusing As atom is 
located at an octahedral interstitial site in ZB GaAs. In 
this interstitial state, the migrating As creates two ad- 
ditional bonds to target tetrahedron apexes, indicating 
that the As atom starts to supply electrons to the target 
tetrahedron. Consequently, the target tetrahedron vol- 
ume further contracts to 0.662Vb. This volume is still a 
bit larger than for a negatively charged As vacancy, but 
notably smaller than for a neutral As vacancy. Passing 
through the AEF face of the target tetrahedron, the dif- 
fusing As atom enters the target vacancy. The diffusion 
process is thus completed. When the As atom diffuses 
into the target tetrahedron, it establishes the As-Ga bond 
to the G-site. The migration from the octahedral intersti- 
tial site to the final state has a shallow migration barrier 
of 0.19 eV which is 0.04 eV lower than the backward dif- 
fusion barrier. We note that the contraction of the AG 
edge of the empty target tetrahedron has moved over to 
the BC edge of the (finally empty) initial tetrahedron. 

Figure [oj^b) shows the V^^ diffusion path in ZB GaAs, 
which is similar to diffusion with a slightly lower mi- 
gration barrier of 1.91 eV. Again, we find the escape bar- 
rier from the initial tetrahedron to be the rate-limiting 
step. In the interstitial space, there is also an interstitial 
As configuration located in a shallow dip, and it occu- 
pies an octahedral interstitial site in ZB GaAs. The as- 
sociated initial and target tetrahedra at this stage have 
volumes of 1.102^0 and 0.673Vb, respectively. We note 
that the tetrahedron volume expands while loosing an 
electron (see Table [l|). The diffusing As atom at the in- 
terstitial site has to overcome an entry barrier of 0.13 eV 
to migrate into the target tetrahedron, which is 0.19 eV 
lower than that of the reverse diffusion process. The ini- 
tial tetrahedron volume is reduced to the V^^ value after 
the diffusing As passes through the second saddle point. 
We note that the two saddle points along the diffusion 
path have different energies. Although an oversimplified 
consideration solely based on the atomic positions sug- 
gests that a symmetric mechanism might exist, the ob- 
served symmetry breaking is physical because the charge 
balance between the two tetrahedra under the presence 
of the migrating As atom is non- symmetric. Therefore, 
it is meaningful to distinguish between the escape (first) 
saddle point or the entry (second) saddle point to be 
rate-limiting. 

Our results qualitatively agree with the previous the- 
oretical study which reported migration barriers of 
2.41 eV and 2.38 eV for V^^ and V^^ diffusion in ZB 
GaAs, respectively, using atomic orbital basis sets. Both 
in their and in our work, the computed diffusion barrier 
of a is slightly lower than that of a V^^. However, 
in our plane- wave-based calculations, diffusion barriers of 
and are computed to be 1.93 eV and 1.91 eV, re- 
spectively, which are about 0.5 eV lower in energy than 
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(a) y+s/wz - ab 



2.5 




FIG. 7: (Color online) The energy path of V^^ (a) and Vl~ 
(b) diffusion in the ab plane of WZ GaAs. The structures of 
two ground state minima (i and v), an interstitial state (in), 
and transition states {ii and iv) are illustrated. The color 
scheme is the same as in Fig. [6] 



the values reported in Ref. 26, We believe that these 
differences are due to the more complete basis set for 
the wave functions in our calculations, which allows for a 
more accurate description of the energetics, in particular 
at transition state geometries. 

Next, we perform analogous calculations for V^^ and 
diffusion in WZ GaAs. Due to the lower crystal sym- 
metry of WZ, one needs to consider two pathways; i.e., 
an As vacancy can diffuse within the ab plane by hopping 
from site 1 to site 2 or along the c axis by hopping from 
site 2 to site 3, as illustrated in Fig.jljb). First, we dis- 
cuss the vacancy diffusion within the ab plane. Figure [7| 
shows that and V^^ can diffuse in the ab plane in 
WZ GaAs with migration barriers of 1.79 eV and 1.60 
eV, respectively. As discussed above, there are two dif- 
ferent minimum configurations for and V^^ defects. 
Ml and M2, in WZ GaAs. The transition from the more 
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TABLE II: Diffusion constants for Vas diffusion in ZB and 
WZ GaAs at T = 700 K. All numbers are given in units of 
cm^/s. The Vas diffusivity along the c direction in WZ GaAs 
is the lowest. 





ZB 


WZ-ab 


WZ-c 


Vas 


2.11 X 10"'^ 
2.94 X 10"^^ 


2.13 X 10"'^ 
4.97 X 10"^^ 


4.46 X lO--"^ 
2.05 X 10"^^ 



stable M2 to the Ml configuration costs relatively low 
activation energies, 0.11 eV for a and 0.17 eV for a 
V^g. We basically find similar diffusion processes as in 
ZB GaAs, i.e., an As escape from the initial tetrahedron, 
an As motion within the interstitial space, and an As 
entry into the target tetrahedron. In the intermediate 
configuration in the ab plane, the As atom is located in 
an octahedral interstitial site in WZ GaAs. Thus, the 
local symmetry of the intermediate is the same for WZ 
and ZB along this path. However, the overall migration 
barriers are up to 0.2-0.3 eV lower in WZ than in ZB. In 
the diffusion, we note that the volume contraction 
of the initial tetrahedron takes place earlier than in the 
other cases (see Fig. [7|^b)), i.e., its volume is reduced to 
0.6911^ when As reaches the interstitial configuration. 
In comparison, for V^^ diffusion in ZB GaAs, the initial 
tetrahedron has a volume of 1.1021^ when the diffusing 
As is at the interstitial site. Moreover, for V^s in WZ 
and ZB crystals, the associated initial tetrahedron vol- 
ume is about 1.8^0 to 1.9Vo when the diffusing As is at 
the interstitial site. This early compression of the initial 
tetrahedron for the V^^ diffusion also results in a com- 
pressed interstitial space. In the interstitial space, the 
diffusing As atom forms three bonds to the target tetra- 
hedron and one bond to the initial tetrahedron. Hence, 
the interstitial is fully coordinated. This could ex- 
plain the low migration barrier of 1.60 eV associated with 
the diffusion within the ab plane in WZ GaAs. 

Secondly, we discuss the Vas diffusion along the c-axis 
in WZ GaAs. Figure [s] shows that and V^g defects 
can diffuse along the c-axis with relatively high migra- 
tion barriers of 2.44 eV and 2.07 eV, respectively. Again, 
the diffusion process consists of the three steps intro- 
duced above: The diffusing As atom escapes from the 
initial tetrahedron cage through the ABC face. It enters 
the neighboring octahedral interstitial area through the 
ABF face and finally hops into the target tetrahedron 
through the AFG face. When the diffusing As passes 
the ABF face, it almost crosses the AF edge. This may 
result in the high migration barriers associated with the 
c-axis diffusion in WZ GaAs. For the diffusion, there 
exists a stable interstitial configuration {Hi) as shown in 
Fig.jsj^b). For this configuration, the diffusing As atom is 
located at a tetrahedral interstitial site. It is fully coor- 
dinated and its energy is only 1.03 eV above the ground 
state i configuration. This may result in the lower V^g 
migration barrier of 2.07 eV compared to the diffu- 
sion as shown in Figjsj 



(a) F+s/WZ - c 




(b) WWZ-c 

2.5 h 



(iv) 




FIG. 8: (Color online) The energy path of (a) and V^^ 
(b) diffusion along the c-axis of in WZ GaAs. The structures 
of various stages {i — v) are illustrated. The color scheme is 
the same as in Fig. [6] 

Finally, we discuss the possibility of reducing the con- 
centration of Vas defects in nanowires by annealing under 
an arsenic atmosphere. We use the migration barrier of 
1.93 eV of a V^^ in ZB GaAs to estimate the required 
annealing time. The diffusivity associated with the Vas 
can be estimated as 

D^aluoexp{-AE/kBT), (3) 

where = 10^^/s is the estimated attempt frequency, /cb 
is the Boltzmann constant, the annealing temperature is 
set to be T = 700 K, AE is the migration energy of an 
arsenic vacancy in GaAs, and ao is the distance between 
two neighboring As atoms in a GaAs crystal and is also 
the hop distance for Vas diffusion. In a nanowire growing 
along the [III] direction, VasS are expected to diffuse to 
the side wall. Using a typical radius of a GaAs nanowire, 
R = 20 nm (Ref. [7j), the annealing time required to re- 
move the VasS can be computed as t = / D = 5.26 
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hours. In WZ GaAs, the reduction of the energy barrier 
by 0.14 eV to 1.79 eV reduces the anneahng time by a 
factor of 10 to 31.3 minutes. Thus, we conclude that the 
VasS introduced in a GaAs nanowire during As-deficient 
growth can be annealed within a reasonable time frame 
at a temperature of 700 K. We hereby summarize the 
diffusivities of and V^^ defects in ZB and WZ GaAs 
crystals for various diffusion processes in Table [ll} 



V. DISCUSSION AND SUMMARY 

Using DFT calculations, we have characterized the 
Vas^ GaAs, and Au defects in GaAs crystals that may 
be introduced in GaAs nanowires grown with the help of 
metal droplets under As-deficient conditions at the inter- 
facial growth zone. The substitutional GaAs, Auas, and 
AuGa defects exhibit similar formation energies and de- 
fect levels in both ZB and WZ GaAs crystals. Moreover, 
a Vas defect behaves as a 'negative /7'-system, switch- 
ing from a in a p-doped material to a V^^ in an 
n-doped material, for both GaAs crystals. In case that 
a Au droplet is used a catalyst for nanowire growth, the 
formation of substitutional Au defects is possible and is 
found to be energetically more favorable than the forma- 
tion of GaAs or Vas defects. Moreover, we have shown 
that it is energetically more favorable by about 1 to 2 eV 
for an Au substitutional defect to replace a lattice Ga 
atom than a lattice As atom in either ZB or WZ GaAs. 
Given that DFT-GGA calculations give a too small band 
gap, the calculated acceptor level of AuQa at Ey +0.22 
eV is found to be in reasonable agreement with the ex- 
perimental value of Ey +0.4 eV in ZB GaAs. 

In Au-free, self-catalyzed growth of nanowires that are 
grown at higher temperatures,!^ the formation of GaAs 
could be expected. In particular, for Ga-rich growth 
conditions and n-type material, the Ga^~ species has 
a low formation energy. This could lead to growth of 
non-stoichiometric, Ga-enriched GaAs.'^'^ The incorpo- 
ration of Ga^~ will counteract the n-doping by delib- 
erately added donor species. This may explain why it 
is rather difficult to obtain n-type conductivity in ZB 
GaAs nanowires, and only high concentrations of SuQa 
have proven successful so far.^'^ Under less Ga-rich growth 
conditions, the formation of UasS could be expected, and 
their diffusion might even play a role for the material 
transport in the growth of a GaAs nanowire. We note 
that As interstitials, which could in principle also con- 
tribute to As mass transport, tend to have a higher for- 
mation energy than VasS, apart from very As-rich condi- 
tions.^^ Since we are interested in the As-deficient growth 
conditions below the nanoparticle, we don't consider the 
As interstitial diffusion here. Our results show that i) 
VasS can diffuse within ZB GaAs with a migration bar- 
rier of about 1.9 eV; and ii) VasS diffuse favorably within 
the ab plane in WZ GaAs with somewhat lower migra- 
tion barriers of 1.6 to 1.8 eV. Based on these results, we 



estimate that it takes about five hours to anneal the VasS 

TABLE III: The Madelung constants (a), the dielectric con- 
stants (e) and the dimension of GaAs supercells (L) for com- 
puting the formation energy as described in Eq. [2] The values 
are listed for both ZB and WZ GaAs supercells. For conve- 
nience, e is given in units of eA~^V~^. The dimensionless 
dielectric constant is obtained by multiplying this value with 
14.4. 





a 


e 


L(A) 


ZB 
WZ 


1.638 
1.641 


1.18 
1.12 


17.29 
12.16 



TABLE IV: The values of the potential off-set corrections 
for Vas (AUvoas), Auas (AVau@As), Auoa (AVAu@Ga), and 
GaAs (AVGa(S)As) defects in WZ and ZB GaAs crystals. All 
numbers are given in units of eV. 





AW@As 


AVauOAs 


AVAuOGa 


AVGa@As 


ZB 
WZ 


-0.05 
-0.09 


-0.01 
-0.02 


-0.04 
-0.02 


-0.05 
-0.09 



at 700 K in a ZB GaAs nanowire, but only 30 minutes in 
a WZ nanowire. Thus, an annealing of GaAs nanowires 
under an arsenic atmosphere could be useful to obtain 
samples with a longer lifetime for the charge carriers. 
However, the diffusivity of VasS in the nanowires is too 
low to contribute substantially to the arsenic supply at 
the nanowire growth zone, in particular in WZ GaAs, 
where the diffusion barrier along the c-axis is higher than 
in ZB GaAs. 
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Appendix: Parameters for computing formation 
energies 

The Madelung constants, the calculated dielectric con- 
stants, and the dimensions of ZB and WZ GaAs crys- 



tals are listed in Table III L represents the distance be- 
tween the nearest neighbor defects within a supercell. It 
is taken to be length of the a-axis of a 216 atom ZB 



and a 96 atom WZ GaAs crystal as listed in Table III 
The dielectric constants are the e^a components of perti- 
nent dielectric tensors. They include both electronic and 
ionic contributions. The values of the potential off- set 
corrections for various defect and impurity calculations 
are listed in Table [TVl 
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